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Abstract. The process of molecule formation in the pri- 
mordial gas is considered in the framework of Friedmann 
cosmological models from redshift z = 10 4 to z = 0. First, 
a comprehensive analysis of 87 gas phase reaction rates 
(both collisional and radiative) relevant in the physical 
environment of the expanding universe is presented and 
critically discussed. On this basis, calculations are carried 
out of the abundance of 21 molecular species as function 
of redshift, for different values of the cosmological param- 
eters f2o, t) and Ho, evaluating consistently the molecular 
heating and cooling due to H2, HD and LiH molecules. 
One of the major improvements of this work is the use of 
a better treatment of H recombination that leads to a re- 
duction of a factor 2-3 in the abundance of electrons and 
H + at freeze-out, with respect to previous studies. The 
lower residual ionization has a negative effect on the chem- 
istry of the primordial gas in which electrons and protons 
act as catalysts in the formation of the first molecules. 
We find that in the standard model (h = 0.67, r)io = 4.5, 
Oo = 1 and [D/H] = 4.3 x 10~ 5 ), the residual fractional 
ionization at z = 1 is [e/H] = 3.02 x 10 -4 , and the 
main molecular species fractional abundances [H2/H] = 
l.lxl0~ 6 , [HD/H 2 ] = l.lxl0~ 3 , [HeH+/H] = 6.2xl0~ 13 , 
[LiH+/H] = 9.4 x 10~ 18 and [LiH/LiH + ] = 7.6 x 10~ 3 . We 
devise a reduced chemical network that reproduces with 
excellent accuracy the numerical results of the complete 
model and allows to follow the chemical compositions and 
the thermal properties of a primordial gas in the presence 
of an external radiation field. Finally, we provide accurate 
cooling functions of H 2 , HD and LiH in a wide range of 
density and temperature that can be conveniently used in 
a variety of cosmological applications. 
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1. Introduction 

The study of molecule formation in the post- 
recombination epoch has grown considerably in recent 
years. Saslaw & Zipoy (1967) and Peebles & Dicke (1968) 
were the first to realize the importance of gas phase 
reactions for the formation of the simplest molecule, H 2 . 
They showed that trace amounts of molecular hydrogen, 
of order 10 -6 -10 -5 , could indeed form via the interme- 
diaries species H+ and H" once the radiation field no 
longer contained a high density of photons with energies 
above the threshold of dissociation (2.64 and 0.75 eV, 
respectively). 

The presence of even a trace abundance of H 2 is of 
direct relevance for the cooling properties of the primor- 
dial gas which, in its absence, would be an extremely poor 
radiator: cooling by Ly-a photons is in fact ineffective at 
temperatures ^ 8000 K, well above the matter and radia- 
tion temperature in the post-recombination era. Since the 
evolution of primordial density fluctuations is controlled 
by the ability of the gas to cool down to low tempera- 
tures, it is very important to obtain a firm picture of the 
chemistry of the dust-free gas mixture, not limited to the 
formation of H 2 , but also to other molecules of potential 
interest. In this regard, Lepp & Shull (1984) and Puy et al. 
(1993) have computed the abundances of H 2 , HD and LiH 
as a function of redshift for various cosmological models. 
Although the final abundances of H 2 and HD agree in the 
two calculations, their evolution with redshift is markedly 
different, since the epoch of formation varies by a factor of 
~ 2. Also, the LiH abundance shows a large discrepancy 
of about two orders of magnitude. 

More recently, Palla et al. (1995) have analyzed the 
effects on the chemistry of the pregalactic gas of a high 
primordial D abundance in the light of the controversial 
results obtained towards high redshift quasars (see e.g. 
Tytler & Buries 1997). They found that the abundance of 
H 2 is rather insensitive to variations in the cosmological 
parameters implied by a factor of ~ 10 enhancement of 
primordial [D /H] , while HD and LiH abundances vary by 
larger amounts. However, the abundance of LiH, obtained 
with simple estimates of the radiative association rate, 
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was largely overestimated. Because of the potential rele- 
vance of the interaction of LiH molecules with the cosmic 
background radiation (CBR) (Maoli et al. 1994), a proper 
treatment of the lithium chemistry was necessary. Dal- 
garno et al. (1996) and Gianturco & Gori Giorgi (1996a, b) 
provided accurate quantum-mechanical calculations of the 
main reaction rates. The chemistry of lithium in the early 
universe has been then studied by Stancil et al. (1996) 
and Bouglcux & Galli (1997). Finally, useful reviews of the 
chemistry of the early universe can be found in Dalgarno 
& Lepp (1987), Black (1991), Shapiro (1992), and Abel 
et al. (1997). The latter two studies, in particular, focus 
on the noncquilibrium H2 chemistry in radiative shocks 
which is thought to be of primary importance during the 
gravitational collapse of density fluctuations (see also An- 
ninos et al. 1997). 

In spite of such a wealth of specific studies, a compre- 
hensive analysis of the subject and a critical discussion 
of the reaction paths and rates are still lacking. To over- 
come this limitation, in this paper we present a complete 
treatment of the evolution of all the molecular and atomic 
species formed in the uniform pregalactic medium at high 
rcdshifts (z < 10 4 ). The structure of the paper is as fol- 
lows: in Sect. 2 we describe the H, D, He, and Li chemistry, 
with a critical discussion of the most important rates; the 
evolutionary models are presented in Sect. 3, and the re- 
sults for the standard model and the dependence on the 
cosmological parameters are given in Sect. 4; Sect. 5 in- 
troduces a minimal model which highlights the dominant 
reactions for the formation of H2, HD, HeH + , LiH and 
LiH+; a comparison with the results of previous studies 
is given in Sect. 6, and the conclusions are summarized in 
Sect. 7. Also, the Appendix provides the collisional exci- 
tation coefficients for HD and H2 and cooling function of 
H2, HD and LiH which are needed for the computation of 
the thermal evolution of the primordial gas. 

2. The chemical network 

The chemical composition of the primordial gas consists 
of e", H, H+, H", D, D+, He, He+, Hc++, Li, Li+, Li", 
H 2 , H+, HD, HD+, HeH+, LiH, LiH+, H+, and H 2 D+. 
The fractional abundaces of these species are calculated as 
function of redshift, starting at z = 10 4 where He, H, D, 
and Li are fully ionized. The network includes 87 reactions 
with the rate coefficients taken from the most recent the- 
oretical and experimental determinations. Reactions rates 
are listed separately for each atomic element in Tables 1 to 
4. In each Table, column 2 gives the reaction, and column 
3 the rate coefficient (in cm 3 s -1 for collisional processes, 
in s _1 for photo- processes) . The gas and radiation tem- 
peratures are indicated by T g and T r , respectively. The 
temperature range of validity of the rate coefficients and 
remarks on the rate are given in column 4, together with 
some remarks on how the rate was obtained. When re- 
action rates were not available, we have integrated the 



cross sections and fitted the results with simple analytical 
expressions (generally accurate to within 10%). The last 
column gives the reference for the rate or the cross section 
adopted. 

The assessment of the accuracy of each rate while in 
principle extremely valuable, is hampered by the large 
number of reactions considered. In this work, we prefer 
to limit the discussion to those rates that are critical for a 
correct estimate of the final molecular abundances, with- 
out omitting any important process. In doing so, we will 
consider the reactions that enter in what will be called 
below the minimal model, i.e. the reduced set of processes 
that are needed to compute the evolution of the main 
molecular species H 2 , , HD, HeH+, LiH and LiH+. In 
particular, for hydrogen we compare our rates to those 
used by Abel et al. (1997) in their extensive compilation. 
For D, He and Li we will refer to rates found in the liter- 
ature as well as in databases. For ease of comparison, we 
show the data from experiments/theory, our numerical fit, 
and the rates adopted by other authors in Fig. 1 for H, 
in Fig. 2 for D and He, and in Fig. 3 for Li. In this way, 
one can immediately appreciate the degree of uncertainty 
associated with each rate. 

2.1. Hydrogen chemistry 

(H1)-(H2) H recombination 

The process of cosmological recombination has been 
analyzed in detail in a number of studies, most recently 
by Jones & Wyse (1985), Krolik (1990), Grachev & 
Dubrovich (1991), and Sasaki & Takahara (1993). The 
improved treatment of stimulated processes during recom- 
bination followed by the last two goups has led to a more 
accurate estimate of the electron abundance at freeze-out. 
Compared to the results of Jones & Wyse (1985) the 
residual ionization degree is a factor 2-3 lower, depend- 
ing on the cosmological model. Of course, the reduced 
number of free electrons and protons has a direct effect 
on the chemical evolution of the gas. In this paper we 
have solved the full equation for the time evolution of the 
electron abundance (Eq. (1) of Sasaki & Takahara 1993) 
using the prescriptions for the rate of radiative recombi- 
nation given by Grachev & Dubrovich (1991). In prac- 
tice, the value for the rate of radiative transition from 
the continuum to the excited states is given by the simple 
power-law fit R c2 = 8.76 x 10- n (l + z)-°- 58 cm 3 s" 1 . The 
inverse rate from the excited levels to the continuum is 
R 2c = 2.41 x 10 15 T r 1 - 5 exp(-39472/T r )i? c2 s" 1 . 
(H3) Radiative attachement of H 

We have fitted the data tabulated by de Jong (1972) 
(triangles), based on the cross section of Ohmura & 
Ohmura (1960). As seen in Fig. 1, our fit (solid line) agrees 
well with the expression used by Abel et al. (1997) (dashed 
line). 

(H4) Radiative detachment of H~ 
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Table 1. Reaction rates for Hydrogen species 



rate (cm 3 s 1 or s x ) 



reaction 



notes 



reference 



HI) 
H2) 

H3) 

H4) 
H5) 

H6) 

H7) 

H8) 

H9) 

H10) 
Hll) 
H12) 
H13) 
H14) 

H15) 

H16) 

H17) 

H18) 

H19) 

H20) 
H21) 
H22) 



H+ + e -> H + 7 
H + 7 -f H++ e 

H + e -> H" + 7 

H" + 7 -> H + e 
H + H -> H 2 + e 



H" + H+ 
H + 
H + H+ 



H+ + c 
-> 2H 



H 



7 



H+ + 7 - H + 



+ H -» H 2 + H 

-> 2H 

2H+ + c 

2 — » H3 + H 

, + +7 
H 2 + H + — ► H 



H 

H+ +e 
H 2 + +7 



H+ + H 



H 2 + e 



H + H" 



H 2 + e 2H + c 
H 2 + 7 -> H+ + e 



H+ +H 
H+ + e - 



- H+ + H 2 
H 2 + H 



H 2 +H+ ^ H+ + 7 
H 3 + + 7 -> H+ + H 



16 2 00 
8823 s ! 
T r 7 



fi c2 
i? 2c 

1.4 x lCr 18 7^ 928 e X p 

1.1 x lu _1 T r 2 - 13 exp(- 

1.5 x 1(T 9 
4.0 x l(T 9 T g - ' 17 
6.9 x 10 _9 T g -°' 35 

9.6 x 10 _7 T g -°' 9 

5.7 x 10- 6 T g -°- 5 + 6.3 x 1(T 8 - 

9.2 x 10 _11 T g - 5 + 4.4 x 10~ 13 T g 
dex[- 19.38 - 1.523 log T g + 
1.118(logT g ) 2 - 0.1269(logT g ) 3 ] 
2.0 x 10 1 T r 159 exp(- 2200°) 
1.63 x I0 7 exp(- 32 f«) r 

6.4 x 10~ 1() 

2.0 x io _7 T- " 



9.0 x lO 1 ^ 148 
2.0 x 10~ 9 
irrelevant 



exp (- 



see text 
see text 

fit 

fit 

T g < 300 
T g > 300, fit 
T g < 8000 
T g > 8000, fit 

fit by PAMS 

1 < T g < 32000, fit 
v = 0, fit 
LTE, fit 

v = 0, fit 
fit 



+ H 3.0 x 10" 



exp 



exp 



2.7 x 10 _8 T g " 1 ' 27 exp 
4.4 x 10~ 1() T"' r " 



2.9 x 10 2 T r 156 exp (-i^ff 
7.7 x 10- 9 cxp (-^) 
4.6 x 10- 6 T g -°' 65 
1.0 x 10~ 16 
irrelevant 




DJ 
DJ 

LDZ 

Po 

MAP 

RP, SBD 
Du 

Ar, St 

KAH 

SDGR 

BO 

TH 

KH 



HMF 

SA 

Co 

OR 

SMT 

Su 

GH 

KH 



Ar: Argyros (1974); BO: Bates & Opik (1968); Co: Corrigan (1965); DJ: de Jong (1972) ; Du: Dunn (1968) ; GH: Gerlich & 
Horning (1992) ; HMF: Holliday et al. (1971) ; KAH: Karpas et al. (1979); KH: Kulander & Heller (1978); LDZ: Launay et al. 
(1991); MAP: Moseley et al. (1970) ; MD: Mitchell & Deveau (1983); OR: ONeil & Reinhardt (1978); PAMS: Peterson et al. 
(1971); Po: Poulaert et al. (1978); RP: Ramaker & Peek (1976); SA: Schulz & Asundi (1967); SBD: Stancil et al. (1993); 
SDGR: Schneider et al. (1994); SMT: Sidhu et al. (1992); St: Stancil (1994) Su: Sundstrom et al. (1994); TH: Theard & 
Huntress (1974) 



The cross section for this reaction is known with great 
accuracy (Wishart 1979). We compare in Fig. 1 our fit 
(solid line) to the the rate tabulated by de Jong (1972) 
(triangles). The expression given by Tegmark et al. (1997) 
is shown by a dashed line. 

(H5) Associative detachment of H and H 

There is substantial agreement between the value of 
this rate coefficient (to within ~ 10%) at temperatures 
below ^10 3 K and the experimental measurement at T = 
300 K performed by Schmeltekopf et al. (1967). At higher 
temperatures earlier results by Browne & Dalgarno (1969) 
indicated a slow rise in the value of the rate, whereas the 
more recent calculations of Launay et al. (1991) show a 



significant decline. An uncertainty of about 50 %, however, 
is still present in the data due to the uncertainty on the 
interaction potential. Our fit (solid line in Fig. 1) is based 
on the results by Launay et al. (1991), and is in reasonable 
agreement with the fit of Shapiro & Kang (1987) to the 
data of Bieniek (1980). On the other hand, Abel et al. 
(1997) favor a rate increasing with temperature for T > 
1000 K, assuming a rate based on the cross section of 
Browne & Dalgarno (1969). 

(H7) Mutual neutralization of H and H + 

This rate is still subject to large uncertainties at low 
temperatures (T <, 2000 K). Our fit in Fig. 1 (solid line) 
is based on the cross sections of Peterson et al. (1971) 
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and Moseley et al. (1970) and the experimental point at 
T = 300 K of Moseley et al. (1970). At temperatures below 
^10 3 K, we have extrapolated as a power-law (oc E~ r ) the 
cross section of Moseley et al. (1970). The rate used by 
Abel et al. (1997) is the same as Dalgarno & Lepp (1987) 
and it is shown by the long-dashed line. Finally the rate 
adopted by Shapiro & Kang (1987), taken from Duley & 
Williams (1984), is shown by the dashed line. 

(H8) Radiative association of H and H + 



This reaction rate has been computed quantum-me- 
chanically by Ramaker & Peek (1976) and more recently 
by Stancil et al. (1993). The two calculations agree to 
within 3% from T = 10 K up to T = 10 6 K. In Fig. 1 we 
show the values tabulated by Ramaker & Peek (1976) (tri- 
angles), our polynomial fit (solid line), and the fit by Abel 
et al. (1997) (dashed line) taken from Shapiro & Kang 
(1987). Since this reaction is important at large redshifts 
when the temperature is well above ^10 3 K, the departure 
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of the fit at low T is of little consequence on the final H2 
abundances. 

(H9) Photodissociation of H^" 

The evolution of Hj, and, to a large extent, that of 
H2 depend crucially on the rate of this reaction. Unfor- 
tunately, calculations of the photodissociation cross sec- 
tion are available only for a sparse set of vibrational levels 
(Dunn f 968, Argyros 1974). Total cross sections, averaged 
over a LTE distribution of level populations, have been 
computed by Argyros (1974) for 2500 K < T < 26000 K 
and by Standi (1994) for 3150 K < T < 25200 K. The 
theoretical calculations are in excellent agreement and re- 
produce satisfactorly well the experimental results by von 
Busch & Dunn (1972). In Fig. 1 we show the LTE pho- 
todissociation rate (solid line), compared with the rates 
for photodissociation from the v = and v = 9 levels 
(dashed and long-dashed lines, respectively). In our stan- 
dard model we adopt the LTE value, but we also show the 
sensitivity of the results to other choices of the reaction 
rate. 

(H10) H2 formation via 

As in most other studies, we assume a rate independent 
of temperature, corresponding to the measurement at T = 
300 K by Karpas et al. (1979). The Langevin rate is shown 
by the long-dashed line. An uncertainty of a factor of <~ 6 
is still present in the literature as indicated by the dashed 
line in Fig. 1 (Culhane & McCray 1995). Such a reduction 
only affects the final Hj abundance by almost the same 
factor, while leaves the H 2 abundance unchanged within 
a few percent. 

(H15) Charge exchange of H2 with H + 

We have integrated the cross section measured by Hol- 
liday et al. (1971) and our fit is shown by the solid line. 
The dashed line shows the rate used by Shapiro & Kang 
(1987), based on a simple detailed balance from the re- 
verse reaction. The rate computed by Abel et al. (1997) 
from the calculations of Janev et al. (1987) shows a dra- 
matic drop at T < 10 4 K (long dashed line in Fig. 1). 
However, Janev et al.'s (1987) data are accurate only for 
T higher than <~ 10 4 K, and should be disregarded at 
lower temperatures. In any case, as noted by Abel et al. 
(1997), large differences in this range of temperatures do 
not affect significantly the final abundances of H 2 and H^ . 

2.2. Deuterium chemistry 

(D1)-(D2) D recombination and ionization 

These processes have been treated in the same way as 
hydrogen (sec (HI) and (H2)). 
(D3)-(D4) Charge exchange of D and H 

These reactions control the evolution of D+ . For (D3) 
we have fitted the tabulated rate coefficient of Watson et 
al. (1978) which shows a weak temperature dependence. 
The rate coefficient for the reverse reaction (D4) is ob- 
tained by multiplying the direct one by exp(43/T). The 
fit is shown in Fig. 2 together with the data by Watson 



et al. (1978) (triangles). Puy et al. (1993) assume a con- 
stant rate that differs from ours by a factor of 2 in the 
temperature interval 10-1000 K. 
(D8) Formation of HD via H 2 

We have used a rate equal to the Langevin collision 
rate constant (solid line in Fig. 2). This rate is slightly 
bigger (2.1 x 10~ 9 vs. 1.7 x 10~ 9 ) than the experimental 
result of Smith et al. (1982) at T = 205 K, but well in 
agreement at T = 300 K (filled circles). These values are 
larger than the previous measurements by Fchsenfcld et 
al. (1973) (filled triangles) and the theoretical estimate by 
Gerlich (1982). 

(D10) Destruction of HD via H+ 

This is the reverse reaction of (D8) which is endoergic 
by 464 K. The solid line in Fig. 2 is the rate given by Smith 
et al. (1982). The dashed line is the theoretical result by 
Gerlich (1982). The data points at T = 205 and 295 K are 
from Henchmann et al. (1981). 

2.3. Helium chemistry 

(He8) Radiative association of He and H+ 

This rate of direct radiative association (He7) is quite 
small compared to inverse rotational predissociation, as 
shown by Flower & Roueff (1979), Roberge & Dalgarno 
(1982) and Kimura et al. (1993). The solid line in Fig. 2 
shows our fit to the photorates computed by Roberge & 
Dalgarno (1982) (triangles). For comparison, the dashed 
line shows the rate of direct radiative association (Roberge 
& Dalgarno 1982). 

(Hell) Destruction of HeH+ by H 

The rate coefficient has been measured by Karpas et al. 
(1979) at T = 300 K (shown by a filled circle in Fig. 2). We 
have used a temperature-independent rate, which is equal 
to about one-half of the Langevin collision rate constant. 
(Hel4) Photodissociation of HeH+ 

This rate has been computed from detailed balance of 
the direct radiative association using a threshold energy 
of 1.96 eV and the rate coefficient of (Hc7). The result 
is shown by the solid line in Fig. 2. For comparison we 
also show the rate for reaction (Hcl5) obtained by fitting 
(dashed line) the results of Roberge & Dalgarno (1982) 
(triangles). 

2.4- Lithium chemistry 

(Li3)-(Li4) Mutual neutralizations of Li and H 

Following Standi et al. (1996), we have adopted the 
same rate for (Li3) and (Li4). The solid line in Fig. 3 
shows the rate obtained using the cross section of Peart 
& Hayton (1994). The long-dashed line is the rate deter- 
mined by Janev & Radulovic (1978). The dashed line is 
the rate adopted by Standi et al. (1996). 
(Li5) Radiative attachment of Li and e~ 

We have used the cross section computed by Ramsbot- 
tom et al. (1994) for the inverse process (Li6) and detailed 
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Table 2. Reaction rates for Deuterium species 



reaction 



rate (cm 3 s 1 or s x ) 



notes 



reference 



Dl) 
D2) 

D3) 

D4) 
D5) 

D6) 

D7) 
D8) 

D9) 

D10) 

Dll) 
D12) 
D13) 
D14) 
D15) 
D16) 
D17) 
D18) 
D19) 
D20) 
D21) 
D22) 
D23) 



D+ + e - 
D + 7 -> 

D + H+ 

D+ + H 



+ D + 7 
D++ e 

-» D + + H 

-> D + H+ 



D + H — » HD + 7 

D + H 2 -> H + HD 

HD+ + H -> H+ + HD 
D+ + H 2 -»■ H+ + HD 

HD + H -> H 2 + D 



HD + H 

HD + H+ 
D + H+ - 
D+ + H - 
HD+ + 7 
HD+ + 7 - 
HD+ + e -> 
HD+ + H 2 
HD+ + H 2 
D + H+ -> 
H 2 D+ + e - 
H 2 D+ + e - 
H 2 D+ + c - 
H 2 D+ + H 2 



-* H 2 + D+ 

H 2 + H 2 D+ 



HD + + 7 
HD+ + 7 

* H + D+ 

♦ H + + D 
. H + D 

H 2 D+ + H 
-> H+ + D 
H 2 D+ + H 
-> H + H + D 
- H 2 + D 
-> HD + H 
H+ + HD 



D24) H 2 D+ + H -f H+ + D 



3.7 x 


10" 


-10^0.28 


3.7 x 


10" 


-10710.28 


1.0 x 


10" 


-25 


9.0 x 


10" 


11 exp ( 


2.1 x 


10" 


-9 



(-* 

76 \ 
« J 



3.2 x 10" 11 exp 

1.0 x 10~ 9 exp 
(2.1 -0.4 log T g ) x 10 



464 



7.2 x 10~ 8 T a 



-1/2 



-1/2 



2.0 x 10 _8 T g -1 
1.0 x 10~ 6 T g " 1/2 x 0.73 
1.0 x 10~ 6 T g " 1/2 x 0.07 
1.0 x 10~ 6 T g 
4.7 x 10~ 9 exp 
5.5 x 10 -10 
2.0 x 10 



x 0.20 



Tg" 1 exp 



(-*) 



see text 
see text 

fit 

fit 

estimate 
T g > 250, fit 
same as YLt + H 



T g > 200 



fit 

same as H + H + 
same as H + H + 
same as H \ + 7 
same as H \ + 7 
fit 

same as Hj + H 2 
same as Hj + H 2 
uncertain 



T E > 10 2 



WCD 

WCD 
LS 

ZM 

KAH 

SAA 

Sh 
SAA 

ASa, MBH 

RP 

RP 



St 
TH 
TH 
ASb 
Da, La 
Da, La 
Da, La 

ASa, MBH, He 
ASb 



ASa: Adams & Smith (1981); ASb: Adams & Smith (1985); Da: Datz et al. (1995); He: Herbst (1982); KAH: Karpas et al. 
(1979) La: Larsson et al. (1996); LS: Lepp & Shull (1984); MBH: Millar et al. (1989); RP: Ramaker & Peek (1976); SAA: 
Smith et al. (1982); Sh: Shavitt (1959); St: Stromholm et al. (1995); TH: Theard & Huntress (1974); WCD: Watson et al. 
(1978); ZM: Zhang & Miller (1989); 



balance arguments to compute the rate shown in Fig. 3 by 
the solid line. The fit by Stancil et al. (1996), corrected 
as in Stancil & Dalgarno (1997), is shown by the dashed 
line. 

(Li6) Photodetachment of Li 

The solid line shows the rate obtained from the pho- 
todetachment cross section evaluated by Ramsbottom et 
al. (1994). 

(Li9) Radiative association of LiH 

This rate has been recently determined quantum-me- 
chanically by Dalgarno et al. (1996) (squares) and Gi- 
anturco & Gori Giorgi (1996a) (triangles). Our fit to the 
results of Dalgarno et al. (1996) is shown by the solid line. 
(LilO) Photodissociation of LiH 

The solid line shows the analytical fit to the rate ob- 
tained by detailed balance of reaction (Li9). 
(Lil8) Radiative association of Li+ and H 



We have fitted (solid line) the results of Dalgarno ct 
al. (1996) (squares). The values obtained by Gianturco & 
Gori Giorgi (1996a) are shown by triangles. The dashed 
line represents the analytic fit by Stancil et al. (1996). 
(Lil9) Radiative association of Li and H + 

As in the case of (Lil8), we have adopted the results 
of Dalgarno et al. (1996) (triangles and solid line). 
(Li25) Photodissociation of LiH + 

The dashed line shows the rate obtained by de- 
tailed balance applied to the radiative association reaction 
(Lil8) adopting the rate by Dalgarno et al. (1996). Our 
fit is shown by the solid line. 

3. Evolutionary Models 

We consider now the chemical and thermal evolution of 
the pregalactic gas in the framework of a Friedmann cos- 
mological model. We start the calculations at z = 10 4 with 
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Fig. 2. Rate coefficients for HD and HcH + chemistry. 

all atomic species fully ionized and we follow their recom- 
bination down to to z = assuming that no rcionizing 
events have taken place in the Universe. 

In order to calculate the abundances of the 21 elements 
included in the network, we have solved the set of coupled 
chemical rate equations of the form 



~TT — k{ OTm TljTlk kd ca t11i + . . . , (1) 



where fcf orm and fcdest are the formation and destruction 
reaction rates discussed in Sect. 2, and rii is the number 
density of the reactant species i. 

The temperature of the CBR is given by T r = Tq(1 + z) 
with T = 2.726 K (Mather et al. 1994). The evolution of 
the gas temperature T g is governed by the equation (see 
e.g. Puy et al. 1993; Palla et al. 1995) 

^ = - 2T £ + ^rd( r - A)co mp to„ + (r - A) mol ]. (2) 
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Table 3. Reaction rates for Helium species 



rate (cm 3 s 1 or s x ) 



reaction 



notes 



reference 



Hel) He ++ + e -» He+ + 7 

He2) He+ + 7 -> He ++ + e 
He3) He+ + e -> He + 7 



He4) He + 7 -» He + + e 
He5) He + H+ -> He+ + H 
He6) He+ + H -> He + H+ 



He7) He + H+ 

He8) He + H 4 

He9) He + H+ 

HelO) He+ + H 



HeH+ + 7 
-> HeH+ + 7 

HeH+ + H 
HeH+ + 7 



Hell) HeH+ + H — > He + 

Hel2) HeH+ + c -> He + H 

Hel3) HeH+ + H 2 -> H+ + He 

Hel4) HeH+ + 7 -> He + H+ 

Hel5) HeH+ + 7 -»■ Hc+ + H 



1.891 x 10" 




1 + 

5.0 x lO 1 ^ 163 exp (- 
3.294 x 10~ n 



15.54 



1 + 



1 + 



3.676 X10 7 



1.0 x 10 4 T r 123 exp (- 



4.0 x W- S 'T* 



-25^2.06 



3.7 x 10~^T & 
1 + 9.9 exp ( 



x 

2570 

O.B 
-1.06 



6717 \ 



5.0 x 10" 
7.6 x 10" 18 T g 
3.45 x 10~ 16 T g ~ 

3.0 x 10~ 10 exp 

1.6 x 10" 14 T g -° 

1.0 x 10~ 15 

9.1 x 10~ 10 

1.7 x 10- 7 T g -°' : 
1.3 x 10~ 9 

6.8 x lO" 1 ^ 16 exp (- 
7.8 x 10 3 T r 12 exp (- 



240000 



15.54 



fit 



fit 

T g > 8000, fit 

fit by KF T g > 6000 

rad. ass. 
T g < 10 3 , 

T g > 10 3 (inv. prediss. 



T g < 4000, fit 
T g > 4000, fit 



fit 
fit 



VF 

Os 

VF 
Os 

KLDD 

ZDKL 
RD, KLDD 

RD 
Bl 

ZD 

KAH 

YM 

Or 

RD 

RD 



Bl: Black (1978); KAH: Karpas et al. (1979); KLDD: Kimura et al. (1993); Or: Orient (1977); Os: Osterbrock (1989); RD: 
Roberge & Dalgarno (1982); VF: Verner & Ferland (1996); YM: Yousif & Mitchell (1989); ZD: Zygelman & Dalgarno (1990); 
ZDKL: Zygelman et al. (1989); 



The first term represents the adiabatic cooling associated 
with the expansion of the Universe, R being the scale fac- 
tor. The other two terms represent respectively the net 
transfer of energy from the CBR to the gas (per unit time 
and unit volume) via Compton scattering of CBR photons 
on electrons 

(r A) 4WT r 4 (T r -T g ) 

(1 - A) Compton — n c , (6) 

and via excitation and de-excitation of molecular transi- 
tions 

(r - A) mo] = ^n k y^^XjCjj - x,( ;,, •//(',,. (4) 

k i>j 

where Cij and Cji are the collisional excitation and de- 
excitation coefficients and Xi are the fractional level pop- 
ulations. For the molecular heating and cooling of the gas 
we have considered the contribution of H 2 , HD and LiH. 



A full discussion of the molecular parameters is given 
in the Appendix, where we present analytical fits and 
plots of the cooling functions in the temperature range 
10K<T g <10 4 K. Given the large range of validity, the 
cooling functions can be used in a variety of cosmological 
applications. 

The energy transfer function (r — A) mo i can become 
an effective heating (cooling) source for the gas if the rate 
of collisional de-excitation of the roto- vibrational molecu- 
lar levels is faster (slower) than their radiative decay (see 
Khersonskii 1986, Puy et al. 1993). In general, however, 
the contribution of molecules to the heating/cooling of the 
gas is very small. 

The chemical network is completed by the equation for 
the rcdshift 



dz H (l + z) 2 y/l + iloz' 
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Table 4. Reaction rates for Lithium species 



rate (cm 3 s 1 or s x ) 



reaction 



notes 



reference 



Lil) Li+ + e -> Li + 7 



Li2) Li + 7 - 
Li3) Li+ + H" 



Li+ + e 

-> Li + H 



Li4) Li" + H + — > Li + H 

Li5) Li + e -> Li + 7 

Li6) LR + 7 -> Li + e 

Li7) Li + H+ -> Li+ + H 

Li8) Li + H+ -> Li+ + H + 7 

Li9) Li( 2 S) + H — > LiH + 7 

LilO) LiH + 7 -> Li( 2 S) + H 

Lill) Li( 2 P) + H -> LiH + 7 

Lil2) Li( 2 P) + H -> LiH + 7 

Lil3) Li + H" -> LiH + e 

Lil4) Li + H -> LiH + e 

Lil5) LiH+ + H -> LiH + H+ 

Lil6) LiH + H+ -> LiH+ + H 

Lil7) LiH + H -> Li + H 

Lil8) Li+ + H -»■ LiH 

Lil9) Li + H+ -»■ LiH 

L120) LiH + H+ -> LiH 

Li21) LiH + H+ -> Li+ + H 2 

Li22) LiH+ + e -> Li + H 

Li23) LiH+ + H -> Li + Hj 

Li24) LiH+ + H -> Li 4 

Li25) LiH+ + 7 -»■ Li+ + H 

Li26) LiH+ + 7 -> Li + H+ 



+ 7 
+ 7 
+ H 



1.036 x 10 _11 [A/Tg/107.7x 
(1 + v / ; T^7l07?7) a612 x 

(1 + ^T g /1.177 X 107^-3881-1 



'+ 



1.3 x l0 3 T r 1 - 45 e X p(-S2ga; 
6.3 x 10" 6 T g -°- 6 - 7.6 x 10 
2.6 x 10" 10 T g - 8 + 2.7 x 10~ 14 T g 

6.1xlO-"2?-«exp(-JL_) 
1.8x 10 2 T r 14 exp(-^) 
2.5xlO-«T»exp(- I | ; ) 

1.7xl0- 13 T- 051 exp (-3^ 
(5.6 x 10 19 T g -° 15 + 
"" 3 exp(- 



— -jt 1 2»2UUU / 

' + 7.2 x lO 15 '/;] 21 ) 1 

8.3 x l0 4 r r °- 3 e X p(-2f20) 

2.0xlO- 16 T°-exp(-^) 

1.9 x 10- 14 T g -°' 34 
4.0 x 10" 10 
4.0 x 10 -10 

1.0 x 10- 11 cxp(-^!2) 

1.0 x 10 -9 
2.0 x 10~ n 
dex[-22.4 + 0.999 log T g - 0.351(log T g ) 21 
4.8 x 10- 14 T g -°' 49 
1.0 x 10" 9 
1.0 x 10 -9 
3.8 x 10- 7 T g -°' 47 

9.0 x l - 10 exp(-^!o'l 



+ H 2 3.0 x 10 _1 ° 



7.0 x 10 2 exp(-i§Pj 
3.35xl0 7 T r exp(- 9 ^) 



det. bal. 
fit 

same as Li + + H~ 

det. bal. 

fit 



det. bal. 

A 1 E + — ^ X J E 4 

s 1 n — ^ X X E+ 

estimate 

estimate 

estimate 

estimate 
estimate 



estimate 
estimate 
estimate 

estimate 

estimate 
det. bal. 
det. bal. 



VF 

PHa 
PHa 

RBB 
RBB 
KDS 
SZ 

DKS, GGGa 

GGGb 

GGGb 

SLD 

SLD 

SLD 

SLD 
SLD 

DKS, GGGa 

DKS 

SLD 

SLD 

SLD 

SLD 

SLD 

DKS, GGGa 
SLD 



DKS: Dalgarno et al. (1996); GGGa: Gianturco & Gori Giorgi (1996a); GGGb: Gianturco & Gori Giorgi (1996b); KDS: 
Kimura et al. (1994); PHa: Peart & Hayton (1994); RBB: Ramsbottom et al. (1994); SLD: Stancil et al. (1996); SZ: Stancil & 
Zygelman (1996); VF: Verner & Ferland (1996) 



where Ho = 100 h km s _1 Mpc -1 is the Hubble constant 
and fio the closure parameter. 

The density n(z) of baryons at redshift z is 

n(z) = n h n cr (l + zf, (6) 

where n cr — 9.19 x 10~ 6 ft. 2 cm -3 is the critical density, 
and fib is the ratio of the total to the critical density. 
This quantity is related to the baryon-to-photon ratio 77 
by n b = 3.66 x lO^io^ 2 , with 7710 = I0 w r). 

4. Results 

4-1. The standard model 

Our standard model is characterized by a Hubble constant 
h = 0.67 (van den Bergh 1989), a closure parameter Oq = 



1, and a baryon-to-photon ratio 7710 = 4.5 (Galli et al. 
1995). The initial fractional abundance of H, D, He and 
Li are taken from the standard big bang nucleosynthesis 
model of Smith et al. (1993). 

The results are shown in Fig. 4, separately for each 
atomic element. The residual ionization fraction at z = 1 
is [e/H] = 3.02 x 10~ 4 . Comparison with Model III of SLD 
which has the same values of the cosmological parameters 
(see Table 5) shows that the freeze-out ionization fraction 
is a factor of 2 smaller in our case, owing to the different 
treatment of H recombination. 

The evolution of the abundance of H2 follows the well 
known behaviour (e.g. Lepp & Shull 1983, Black 1991) 
where the initial steep rise is determined by the chan- 
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Fig. 3. Rate coefficients for LiH chemistry. 



nel (reactions H9 and H10) followed by a small contribu- 
tion from H~ at z ~ 100 (reactions H3, H4 and H5). The 
freeze-out value of H2 is [H2/H] = 1.1 xlO -6 , a value some- 
what lower than found in similar studies, as a consequence 
of the net reduction of the ionization fraction. 

Fig. 4 shows that all the ionic species reach very low 
fractional abundances. In particular, the steady drop of 
H at z < 100 is determined by mutual neutralization 
with H+ (reaction H7). As discussed before, the rate for 
this reaction is still uncertain, but the net effect on the 



final abundance of H2 is quite small. The reason is that 
the rates differ only at low temperatures (below ~ 100 
K) when the H2 formation via the H channel is already 
completed. In fact, adopting the Dalgarno & Lepp (1987) 
rate for reaction (H7) leads to a 10% difference in the 
H2 abundance. On the contrary, the abundance of H~ is 
quite affected by the choice of the rate of (H7): our final 
H~ abundance at z = 1 is three times smaller than that 
obtained using Dalgarno & Lepp (1987). 
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redshift (1+z) 



Fig. 4. The evolution of all the chemical species considered in the standard model as a function of redshift. The four 
panels show the results for H (upper left), D (upper right), He (lower left), and Li (lower right). 



The formation of Hj is controlled by radiative associa- 
tion (H8) and photodissociation (H9) at z > 100, whereas 
at lower redshifts there is a contribution from HeH + (re- 
action Hell) that explains the gentle rise up to [HJ/H2] 
~ 10 -6 . An additional channel for H^ formation is the as- 
sociative ionization reaction H+H(n = 2) — ► H^+e (Rawl- 
ings et al. 1993; see also Latter & Black 1991 for H 2 for- 
mation). The rate coefficient for this process is ~ 4 orders 
of magnitude faster than (H8) but the extremely low frac- 



tional population of the n = 2 level of H (<C 10~ 13 ) makes 
associative ionization not competitive with radiative asso- 
ciation. The evolution of H^" and, to a large extent, the 
abundance of H2 depend crucially on the adopted pho- 
todissociation rate. This is clearly shown in Fig. 5 where 
the results obtained with various choices of the rate dis- 
cussed in Sect. 2.1 are compared with the results of the 
standard model (shown by the solid line). Photodissocia- 
tion of H^ from v = (dashed line) would result in an 
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Table 5. Ionization fraction at z = 10 



#0 


O 




T 


[e/H] 


[e/H] 








(K) 


(this work) 


(other) 


67 


1 


0.0367 


2.726 


3.3(-4) 


6.5(-4) [SLD-III] 


50 


0.1 


0.1 


2.7 


5.1(-5) 


5.0(-5) [B] 


50 


1 


0.1 


2.726 


1.7(-4) 


3.2(-4) [SLD-I] 


100 


1 


0.1 


2.7 


8.3(-5) 


7.8(-5) [LB] 



B: Black (1991); LB: Latter & Black (1991); SLD-I, SLD-III: 
Models I and III of Standi et al. (1996) 



enhancement of a factor <~ 200 of the final H2 abundance, 
whereas photodissociation from v = 9 (long-dashed line), 
would delay the redshift of formation of Hj from z ~ 10 3 
toz~ 300 and of H 2 from z ~ 600 to z ~ 250. How- 
ever, the asymptotic abundances of both species are left 
unchanged. 

We believe that the uncertainty in the redshift evolu- 
tion of H 2 is not as large as that shown in Fig. 5, because 
photodissociation from v = is extremely unlikely to oc- 
cur at z ~ 10 3 . The actual behaviour of the H 2 abun- 
dance should be intermediate between the two curves in 
Fig. 5 corresponding to photodissociation from v = 9 and 
LTE. The conditions for a distribution of level populations 
in LTE may be marginally satisfied, considering the fact 
that the radiative decay rates (Posen et al. 1983) are much 
faster than the collisional excitation rates at the typical 
densities, but we know from the work of Ramakcr & Peek 
(1976) that the reverse reaction (H8) forms Hj preferen- 
tially in excited states. In such a case, photodissociation 
would take place from high vibrational levels and the LTE 
rate may underestimate somehow the actual rate. Thus, 
we conclude that the uncertainty is limited to the red- 
shift of formation of H 2 molecules and not to their final 
abundance. 

As for deuterium, Fig. 4 shows that the evolution of D + 
is very sensitive to the charge exchange reactions (D3)- 
(D4). At z < 40, the abundance of D+ drops precipi- 
tously due to the exponential factor in the reaction rate 
(D3). The only molecule formed in significant amount is 
HD, whose evolution with redshift follows closely that 
of H 2 . However, its abundance at freeze-out is [HD/H 2 ] 
~ 10~ 3 . More important, the asymptotic value of [HD/H 2 ] 
becomes about 10~ 3 with an enhancement factor of ~30 
from the initial [D/H] abundance. This is due to the large 
fractionation implied by reaction (D8). The behaviour of 
HD+ follows that of Hj up to z ~ 100, while at lower red- 
shifts its evolution is dominated by reaction (D12). The 
flattening at z < 10 reflects the behaviour of the reaction 
rate at temperatures below ~ 20 K. 

The case of H 2 D+is of interest after the suggestion by 
Dubrovich (1993) and Dubrovich & Lipovka (1995) that 



chemistry of the early Universe 



,,,,,,, 






H 






H + 




/ / " 
/ / 
/ / 


/ 




/■'' ft 






/ / / 

n / / 7\ 


hj ; 




< / ' 
,7 / / 

III 





4 3 2 1 

log (1 + z) 

Fig. 5. Effects of varying the photodissociation rate of H^ 
on the evolution of H 2 and H^ . The solid curves are for 
the standard models, whereas the results obtained with 
photodissociation from the v = and v = 9 vibrational 
levels of H^ are shown by the dashed and long dashed lines, 
respectively. 

in the early universe this molecule could reach significant 
abundances (up to [H 2 D + /H 2 ] = 10 -5 ) and that spec- 
tral features associated with the rovibrational transitions 
might be detectable in the spectrum of the CBR. Such a 
high abundance would result from the complete conver- 
sion of H^" into H 2 D + because of deuterium fractionation 
at low temperatures (T < 50 K). However, it is very un- 
likely that this mechanism could work in the conditions of 
the primordial gas, since the abundance of H3 ions is quite 
low due to the absence of ionizing sources (see Fig. 4) . In 
present-day molecular clouds, searches for the rotational 
transitions of H 2 D + have yielded a fractional abundance 
of 3 x lO" 11 (Boreiko & Betz 1993). The results shown 
in Fig. 4 indicate that the final abundance of H 2 D + is ex- 
tremely small [H 2 D+/H] ~ 5 x 10~ 14 , and [H 2 D+/H+] 
~ 10~ 2 . It is clear from our models that the presence of 
H 2 D+ molecules can hardly affect the shape of the CBR 
(in order to see detectable effects on the CBR, Dubrovich 
& Lipovka (1993) assume an abundance of 10~ 8 ). On the 
other hand, the isotopic ratio is greatly enhanced with re- 
spect to the primordial abundance [D/H] = 4 x 10~ 5 , due 
to fractionation effects. 

According to the results shown in Fig. 4, the main 
molecular species containing helium is HcH + , formed by 
the radiative association of He and H + . As emphasized 
by Lepp & Shull (1984), this reaction is slower than the 
usual formation mode via association of He + and H, but 
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since the abundance of He + is quite small, reaction (Hc8) 
takes over. The HeH + ions are removed by CBR photons 
and at z ^ 250 also by collisions with H atoms to reform 
(reaction Hell). This explains the abrupt change in 
the slope of the curve of HeH + shown in Fig. 4. The hnal 
abundance is small, [HeH+/H] ~6x 10~ 13 at z = 1. 

Finally, the chemistry of lithium is complicated, but 
the molecular abundances arc indeed very small. The more 
abundant complex is LiH+ whose formation is controlled 
by the radiative association of Li + and H (see Dalgarno 
& Lepp 1987). As shown in Fig. 4, lithium remains more 
ionized at low redshifts and this explains why LiH is less 
abundant than LiH + . The final abundance of LiH+ is ~ 
1(T 17 at z = 1. 

4-2. Dependence on the cosmological parameters 

The final values of the molecular abundances depend on 
the choice of the cosmological parameters. Each model 
is in fact specified by three parameters 7/10, flo and h. 
The observational (and theoretical) uncertainties associ- 
ated with each of them are still rather large. We have 
selected as standard model that that better fits the con- 
straints imposed by the observations of the abundances 
of the light elements. To gauge the effects that a change 
of f?io induces on the chemistry network, we have run a 
model with 7710 =8.0, the maximum value still compatible 
with the constraints on standard big bang nucleosynthesis 
given by observations of lithium in Pop-II stars (Bonifa- 
cio & Molaro 1997). The results are given in Table 6 for a 
direct comparison with the standard case. Note that the 
abundance of H2 is the same in the two models, making 
this molecule a poor diagnostic of cosmological models. 
Lithium molecules are quite sensitive to 7710, showing an 
increase of a factor of 3-4. For the other species, however, a 
higher value of -qio implies a lower final abundance as a re- 
sult of the lower initial value of each element predicted by 
the standard big bang nucleosynthesis which compensates 
for the higher total baryon density. As shown by Palla et 
al. (1995), larger variations of the molecular abundances 
than those shown in Table 6 can be expected for more 
drastic variations of the other cosmological parameters: 
both LiH and HD can vary by 2 or 3 orders of magnitudes 
by allowing changes to fio an d h of a factor of 5 and 2, 
respectively. However, current observations indicate that 
a high value of the Hubble constant is unlikely, and that 
a universe with f2 < 1, as in the case of a nonzero value 
of the cosmological constant, is less favored by theoretical 
models. 

5. The minimal model 

The chemical network discussed in Sect. 2 consists of 87 
reactions, 69 collisional and 18 radiative. Although the 
integration of such system does not present particular dif- 
ficulties (apart for the intrinsic stiffness of the rate equa- 



tions) or require exceedingly long computer times, for a 
better understanding of the chemistry of the primordial 
gas it is more convenient to reduce the reactions to only 
those that are essential to accurately model the forma- 
tion/destruction of each molecular species. Such a reduced 
system has been called the minimal model. Abel et al. 
(1997) have devised a reduced network in their discussion 
of the non-equilibrium effects and chemical dynamics of 
the primordial gas. However, their application differs from 
ours, since they are mainly interested in the formation of 
H2 in the post-shock layer of a cosmological pancake and 
consider only collisional processes. 

Our minimum model consists of the 33 reactions indi- 
cated in the diagram shown in Fig. 6. The main reactions 
are listed in bold face in Tables 1 to 4. Notice that in or- 
der to describe fully the H 2 chemistry only 11 reactions 
out of 22 arc needed. The formation of H2 involves two re- 
action sequences with H~ and H^ ions whose abundance 
is restricted by photodestruction processes. However, for 
H~ a competing destruction channel is the mutual neu- 
tralization with H + ions which contributes ~ 20-40% to 
the total rate at redshifts z <, 100. The H 2 system would 
be closed if not for the external source of H^ ions coming 
from the destruction of HeH + by H atoms (reaction Hell) 
at redshifts z ^> 250. It is important to consider this path 
since it modifies substantially the evolution of H^~ (see 
Fig. 4) , although it does not affect the final abundance of 
H 2 molecules. As we have already noted, in the absence 
of this channel the final abundance of H^~ would be of the 
same order as that of H3 , because of the importance of 
reaction (H19). 

The chemistry of deuterium is rather simple and only 
6 reactions out of 24 need to be considered. The molecule 
HD is formed rapidly by reaction (D8) which involves ion- 
ized deuterium. Thus, it is important to consider all the 
reactions that determine the ionization balance. The main 
destruction route is by reaction (D10), while other reac- 
tions involving H;j~ (Dll) and H (D9) are much less im- 
portant. 

Finally, the chemistry of Li is more complex than the 
rest and even the reduced network contains a large number 
of reactions (14 out of 26). This is mainly due to the fact 
that lithium remains partly ionized and the routes to the 
formation of LiH and LiH + compete effectively with each 
other. Note also that at redshifts z ^ 100, the main route 
to LiH formation is from associative detachment of Li~ 
ions with H. This process gradually removes Li - and the 
abundance of LiH reaches a constant value. 

6. Comparison with previous work 

In order to make accurate comparisons with the results 
obtained in previous studies, we have run models with 
the same reaction rates adopted in our standard case, but 
differing cosmological parameters to reproduce exactly the 
choices of the various authors listed in Table 7. 
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Table 6. Sensitivity of abundances to -q w 
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0.926 
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2.3(-19) 



Table 7. Values of the cosmological parameters 
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Fig. 6. The minimal model for the main chemical species 

The results are shown in Fig. 7, where we compare the 
predictions of our model (solid lines) with those of Palla 
et al. (1995) (upper left panel), Lepp & Shull (1984) (up- 
per right panel), Puy et al. (1993) (lower left panel), and 
Black (1991) (lower right panel), all indicated by dashed 
lines. A common feature in the first three cases is a dra- 
matic reduction of the abundance of LiH by 7-9 orders 
of magnitude, due to the replacement of the semiclassi- 
cal estimate of the rate of radiative association (Lepp & 
Shull 1984) with the quantal calculations by Dalgarno et 
al. (1996) and Gianturco & Gori Giorgi (1996a). 

With respect to our earlier study of primordial chem- 
istry, the more accurate treatment of H recombination re- 
sults in a decreased residual electron fraction by a factor 
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PGS: Palla et al. (1995); LS: Lepp & Shull (1983); PALLP: 
Puy et al. (1993); B: Black (1991); GS: Giroux & Shapiro 
(1996) 



2-3 (see also Table 5). This, in turn, has a direct impact 
on the H 2 and HD abundances, which are lower by a fac- 
tor <~ 5. Otherwise, the evolution with redshift of these 
two species follows the same behaviour. 

Similar considerations apply to the comparison with 
the results obtained earlier by Lepp & Shull (1984) for H, 
H + and H2. However, it is worth noticing that in their 
model the formation of HD occurs through the analogues 
of the H~ and Hj channels, at rates diminished by the cos- 
mological D/H ratio, plus a negligible contribution from 
direct radiative association of H and D. In contrast, we 
have found that the main reaction responsible for the for- 
mation of HD is the isotope exchange reaction (D8) which 
has no counterpart in the hydrogen chemical network. 
Therefore, despite their higher residual electron fraction, 
their asymptotic abundance of HD is underestimated by 
about one order of magnitude. 

Puy et al. (1993) obtained abundances of H 2 and HD 
about one order of magnitude larger than ours, owing to 
their unrealistic choice of the photodissociation rate (H9) 
by von Bush & Dunn (1968), valid for a LTE distribu- 
tion of level populations at T = 100 C. The higher abun- 
dance of Hj results in a higher abundance of H 2 formed 
through this channel, and, in turn, to a correspondingly 
larger abundance of HD through reaction (D8). 

Finally, the results of Black (1991), Latter & Black 
(1991), and Giroux & Shapiro (1996) allow a direct com- 
parison of the evolution of some less abundant species 
like H~, H^ and HcH + . The agreement with our results 
and those by Black (1991) is extremely good: the abun- 
dance of H, H+ and H 2 agree within less than a factor 
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1.5 from z = 10 3 to z = 10, while the abundances of H~, 
Hj, although differing by a factor ~ 3, follow remark- 
ably the same behaviour with redshift. Only HcH+ shows 
a marked (but unconsequential) difference for z > 300, 
probably due to a different choice of the photodissocia- 
tion rate (Hcl4). On the other hand, the comparison with 
Giroux & Shapiro (1996) shows a different behaviour of 
at z ^ 100 resulting from the neglect of the HeH + 
channel for H^ formation. 



7. Conclusions 

The main results of the present study can be summarised 
as follows: 

1) We have followed the chemical evolution of the pri- 
mordial gas after recombination by computing the abun- 
dances of 21 species, 12 atomic and 9 molecular, by using 
a complete set of reaction rates for collisional and radia- 
tive processes. The rates which are critical for a correct 
estimate of the final molecular abundances have been anal- 
ysed and compared in detail. 

2) One of the major improvements of this work is the 
use of a better treatment of H recombination that leads to 
a reduction of a factor 2-3 in the abundance of electrons 
and H + at freeze-out, with respect to previous studies. The 
lower residual ionization has a negative effect on the chem- 
istry of the primordial gas in which electrons and protons 
act as catalysts in the formation of the first molecules. 

3) In the standard model (h — 0.67, 7710 = 4.5, Clo = 1 
and [D/H] = 4.3 x 10~ 5 ), the residual fractional ioniza- 
tion at z = 1 is [e/H] = 3.02 x 10~ 4 , and the main 
molecular species have fractional abundances [H 2 /H] = 
1.1 xl0~ 6 , [HD/H 2 ] = l.lxlO" 3 , [HcH+/H] = 6.2xl0~ 13 , 
[LiH+/H] = 9.4 x 10~ 18 and [LiH/LiH + ] = 7.6 x 10~ 3 . 

4) As for molecular hydrogen, its final abundance 
does not depend on the model parameters, making this 
molecule a poor diagnostic of cosmological scenarios. The 
largest uncertainty resides in the accurate knowlwedge of 
the photodissociation rate of H^. A detailed treatment of 
the reaction kinetics of this reaction would be required. 

5) We have presented a minimal model consisting of 
11 reactions for H 2 , 6 for HD, 3 for HcH+ and 14 for LiH 
which reproduces with excellent accuracy the results of 
the full chemical network, regardless of the choice of the 
cosmological parameters. 

6) Finally, we have computed accurate expressions for 
the cooling functions of H 2 , HD and LiH in a wide range 
of density and temperature that can be conveniently used 
in a variety of comological applications. 
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Appendix A: Appendix: Excitation Coefficients 
and Cooling 

We have evaluated the contribution of H2, HD, and LiH to 
the heating/cooling properties of the primordial gas. For each 
molecule, the steady-state level populations are obtained at 
any z by solving the balance equations 



,ji (Rj'j + Cj'j 



(A.l) 



where j and j' indicate a generic couple of rovibrational levels. 
The collisional transition probabilities Cjji and Cjij are ob- 
tained by multiplying the corresponding excitation coefficients 
jjji and jj/j by the gas density. The terms Rjji and Rjij are 
the radiative excitation and de-excitation rates, that can be 
expressed in terms of the Einstein coefficients A,y and B 



A n' +Bjj'it(^y,Tr), f <j, 
B, y u(z/,y,T r ), f>3, 



n' 1 



(A.2) 



where uivjji , T r ) is the energy density of the radiation per unit 
frequency at the temperature T r , 



8Trhv : -,, 

u(v jr ,T r ) = j^- [exp(his jf /kT T ) - I]" 1 



(A.3) 



To evaluate the cooling function A mo i(z), the radiative 
and collisional transition probabilities and the frequencies of 
the roto-vibrational transitions must be specified for each 
molecule. These data can be easily found in the literature only 
for H2, that is obviously the best studied case, but are quite 
sparse for the two other species we have considered. As for 
LiH, the relevant quantities have been collected and critically 
examined by Bougleux & Galli (1997), and we refer the inter- 
ested reader to their paper (in particular, see their Appendix 
B). Here, we limit ourselves to a discussion of the adopted 
coefficients for HD and H2. 



A.l. H 2 

Ortho- and para-H2 were assumed to be in their equilibrium 
(3:1) ratio. The radiative rates were calculated by Turner et 
al. (1977). Since in some cases where H2 plays a major role the 
gas temperature can reach more than one thousand degrees, 
it is not sufficient to consider only the ground (v = 0) vibra- 
tional state. The choice of the rovibrational H-H2 collisional 
coefficients is very delicate, since the published data may vary 
up to factors of 50-100 according to different authors (see e.g. 
Shull & Beckwith 1982, Mandy & Martin 1993, Sun & Dal- 
garno 1994, Flower 1997a,b), and the resulting H-H2 cooling 
functions may differ by considerable amounts. For cosmologi- 
cal problems, one generally adopts the analytical expressions 
given by Hollenbach & McKee (1979, 1989) or Lepp & Shull 
(1983). These formulae appear now to suffer considerably from 
the limitations and uncertainties associated with the collisional 
coefficients available at the time. At low temperatures, the un- 
certainty is associated with the choice of the interaction po- 
tential, since it is very difficult, even nowadays, to calculate 
the potential in the interaction region to the requisite level 
of accuracy (~ 10~ 3 hartrees). The recently developed poten- 
tials LSTH (Truhlar & Horowitz 1978), DMBE (Varandas et al. 
1987), BKMP1 (Boothroyd et al. 1991), BKMP2 (Boothroyd et 



16 



Danicle Galli and Francesco Palla: The chemistry of the early Universe 




Fig. 7. Comparison between the results of the standard model (solid lines) and those of Palla et al. (1995) (upper 
left), Lcpp & Shull (1984) (upper right), Puy et al. (1993) (middle left), Black (1991) (middle right), and Giroux & 
Shapiro (1996) (lower left), all indicated by dashed lines. 



al. 1996) and PBSL (Partridge et al. 1993) have been adopted 
by various groups to generate sets of collisional cross sections 
and rate coefficients, using different methods. Table Al shows 
a comparison of the resulting values of 720 at T g = 10 3 K. 

At temperatures higher than ~ 600 K, we have adopted the 
set of rovibrational H-H2 rate coefficients computed by Mandy 
& Martin (1993). The corresponding cooling function has been 
determined by Martin et al. (1996). At lower temperatures, 
we have adopted the collisional coefficients recently computed 



by Forrey et al. (1997) that match those of Mandy & Martin 
(1993) at T = 10 3 K. We have fitted the tabulated results by 
expressions like 

U, j' = a + aiT g + a 2 T g 2 + a 3 T g 3 , (A.4) 

where 7j,j' is in cm 3 s _1 , and the coefficients a, are listed in 
Table A2. Different choices (like, e. g. Flower 1997a,b) are of 
course possible, and produce differences of a factor ~ 2 in the 
cooling function. 
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Table A.l. H-H 2 Collisional Deexcitation: Com- 
parison of Different Results at T„ = 10 3 K 



Potential 
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720 




Reference 






(cm 3 s- 1 
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GT: Green & Thrular (1979); MM: Mandy & Martin (1993); 
LBD: Lepp et al. (1995); SD: Sun & Dalgarno (1994); Fa: 
Flower (1997a); Fb: Flower (1997b); FBDL: Forrey et al. 
(1997) 

Table A. 2. H-H 2 Collisional De-excitation Coef- 
ficients (from Forrey et al. 1997) 
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1.21(-15) 


2.16(-19) 
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Fig. A.l. Cooling function of H2 for H-H2 collisions for 
two values of the density. The thick solid curves show 
the results of our fit as given by Eq. A5-A7. The other 
curves give the cooling functions computed by Lepp & 
Shull (1983) (LS: dot-dash), Hollenbach & McKee (1989) 
(HMK: dash), Martin et al. (1996) (MSM: dot- long dash), 
and Forrey et al. (1997) (FBDL) (thin solid). 



It is convenient to express the H2 cooling function Ah 2 (in 
erg cm 3 s"" 1 ) in the form given by Hollenbach & McKee (1979), 



A H >(H),T g ] 



A(LTE) 



1 + 



7n(H)] ' 



(A.5) 



where Ah 2 (LTE) is the LTE cooling function given by Hollen- 
bach & McKee (1979), n cr is the critical density defined as 



A H2 (LTE) 



n(H) A Hj [n(H)-»0]' 



(A.6) 



and A[n(H) — > 0] is the low-density limit of the cooling func- 
tion, independent of the H density. The latter can be computed 
from the collisional and radiative deexcitation coefficients de- 
scribed above, and the result is well approximated by the ex- 
pression 

logA H2 KH) -> 0] = -103.0 + 97.59 log T g - 48.05(logT g ) 2 
+ 10.80(logT g ) 3 - 0.9032(logT g ) 4 , (A.7) 

over the range 10 K < T g < 10 4 K. 

The H-H2 cooling function is shown in Fig. Al for two val- 
ues of the H density, n(H) = 0.1 cm" 3 and n(H) = 10 6 cm" 3 . 
The cooling functions given by Hollenbach & McKee (1989) 
and Lepp & Shull (1983), also shown in Fig. Al, differ sig- 
nificantly from that computed in this work especially in the 
low-density limit, at both low and high temperatures (see dis- 
cussion in Martin et al. 1996). 

In considering the cooling by H2 it is important not to 
neglect the effect of reactive collisions of H2 with H or H + 



that cause ortho-para interchange and lead to a ratio out of 
the equilibrium value 3:1. Calculations of the relative colli- 
sional coefficients have been made for H2-H + by Dalgarno et 
al. (1973) and Flower & Watts (1984), and for H 2 -H by Sun 
& Dalgarno (1994). In the conditions of the primordial gas 
with H+/H £ 10~ 4 the ortho -para exchange is dominated by 
collisions with protons, and cooling at low temperatures may 
occur primarily through the J = 2 — > transition of para- 
H2. A complete calculation of the cooling function should take 
into account this effect that has been neglected in the results 
presented in this section. 

A.2. ED 

Due to its small dipole moment, D — 8.3x 10 -4 debyes (Abgrall 
et al. 1982), the energy spacing of the rotational levels of HD 
is quite large. The energy levels and the rotational radiative 
transition probabilities have been computed by Dalgarno & 
Wright (1972) from the experimental data of Trefler and Gush 
(1968). The calculations have been extended by Abgrall et al. 
(1982), who considered both dipole and quadrupole transitions 
and included a large number of rotational and vibrational lev- 
els. The two sets of results differ by a factor of two, and we 
have used the most recent determination. 

The coefficients for inelastic scattering of He-HD have been 
computed by Green (1974) using the Shafer & Gordon (1973) 
potential of H2-He, and more recently by Schaefer (1990) 
with the Schaefer & Kohler (1989) potential at temperatures 
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Table A. 3. HD-He Collisional De-excitation Co- 
efficients (from Schaefer 1990) 



10" 19 e 
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8.7 x 10" 


-14 


1.03 
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3.2 x 10" 


13 


1.3 x 10" 


-15 


1.47 



T g < 600 K for < J < 3 and A J = +1,-1-2. De-excitation 
coefficients depend approximately linearly with temperature. 
We have fitted the numerical data of Schaefer (1990) with ex- 
pressions of the form Jjij = ao + Oi Tg , and fit coefficients are 
tabulated in Table A3. These values should be divided by 1.27 
to account for the difference between the (HD, H) and (HD, 
He) interaction potentials (Wright & Morton 1979). 

A. 3. A total cooling function for the primordial gas 

In Fig. A2 we summarize the different contributions to the 
cooling of a gas of primordial composition. Although the curves 
have been obtained in the limits n(H) — + 0, they are valid for 
n(H) Ss 10 2 cm -3 . At larger densities, the cooling functions 
rapidly approach the corresp ond ing LTE values; at interme- 
diate densities, the formula (A.5) represents a very good ap- 
proximation. The cooling functions for H-H2 and H-HD are 
obtained as described in this work, the H-LiH cooling function 
is computed H^" cooling function for collisions with H and e~ 



was computed by Suchkov & Shchekinov (1978). 

The cooling function of HD in the low-density limit for 
T g Si 10 3 K is accurately reproduced by the expression 



Ahd [71(H) — > 0] 



= 2 7 io£ioexp(-£VfcT g ) 
+ (5/3)721^21 cxp(-£ 21 /fcT g 



(A.8) 



where E\o = 128fc, E21 = 255fc, and the 7jj/ are given 
in Table A3. As for LiH, the cooling function in the low- 
density limit is well approximated by the polynomial expres- 
sion logA LiH (nH) -» 0) = c + ci(logTg) + c 2 (logT g ) 2 + 
c 3 (logT g ) 3 + c 4 (logTg) 4 in the temperature range 10-10 3 K. 
With ALiH in erg cm 3 s -1 , the fit coefficients are Co = —31.47, 
ci = 8.817, c 2 = -4.144, c 3 = 0.8292, c 4 = -0.04996. 

In the standard model, the radiation temperature is higher 
than the matter temperature and molecules are a net heating 
source for the gas. Molecular hydrogen dominates the heat- 
ing of the gas at temperatures higher than ~ 150 K, whereas 
HD dominates at lower temperatures. The ability of HD to 
heat/cool a gas when H 2 molecules become unefficient is a re- 
markable property of a gas of primordial composition, and has 
been verified in a number of studies (Varshalovich & Kher- 
sonskii 1977, Shchekinov 1986, Bougleux & Galli 1997, Puy & 
Signore 1997). 
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